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Abstract 

We present a unified treatment of explicit in time, two level, second 
order resolution, total variation diminishing, approximations to scalar 
conservation laws. The schemes are assumed only to have conservation form and 
incremental form. We introduce a modified flux and a viscosity coefficient and 
obtain results in terms of the latter. The existence of a cell entropy 
inequality is discussed and such an equality for all entropies is shown to 
imply that the scheme is an E scheme on monotone (actually more general) 
data, hence at most only first order accurate in general. Convergence for 
TVD-SOR schemes approximating convex or concave conservation laws is shown by 
enforcing a single discrete entropy inequality. 
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Introduction 


Recently there has been an enormous amount of activity related to the 
construction and analysis of finite difference approximations which 
approximate nonlinear hyperbolic systems of conservation laws and which are 
supposed to have the following properties: 

(1) Limit solutions which satisfy a geometric and/or analytic entropy 
condition. 

(2) A bound on the variation of the approximate solutions at least in the 
scalar and linear systems case. This bound is such as to imply the 
absence of spurious oscillations in the approximate solutions. 

(3) At least second-order accuracy in regions of smoothness, except for 
certain isolated points as described below. 

Some examples of the successful computational consequences of this 
activity can be found in the proceedings of the sixth AIAA Computational Fluid 
Dynamics Conference, and elsewhere see e.g., the bibliography in [21], 

Some of the earliest work in the design of schemes having properties (2) 
and (3) above was done by Van Leer [27], [28], There he introduced the 

concepts of flux limiters and higher order Riemann solvers. Recently Harten 
[10], [11] obtained conditions which he showed to be compatible with second- 

order accuracy, and which guarantee that a scalar one-dimensional scheme is 
TVD—total variation diminishing. He constructed a scheme having that 
property, and formally extended it to systems using a field-by-field limiter 
and Roe's decomposition [22], 

We would also like to mention the work of Boris and Book [1] concerning 
FCT schemes. They also used flux limiters to suppress oscillations in their 


schemes 
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Harten's construction in [10] was done first for a fully discrete, 
explicit in time approximation. P. Sweby [24] has investigated the properties 
of various limiters in this context. 

We shall use the term "high resolution scheme" to mean a formal extension 
to systems via a field-by-field decomposition of a scalar, higher than first- 
order accurate, variation diminishing scheme. These schemes do not, in 
general satisfy the entropy condition e.g. expansion shocks exist as stable 
solutions of high resolution schemes based on Roe's (unmodified) scheme. In 
[20], Osher's decomposition and certain flux limiters were used to prove that 
limit solutions to a class of semi-discrete, time continuous high resolution 
schemes do satisfy the entropy condition for hyperbolic systems of 
conservation laws. Convergence of other classes of semi-discrete, time 
continuous, high resolution approximations to scalar convex conservation laws 
was proven in [19] and [20]. 

A systematic recipe for constructing serai-discrete high resolution schemes 
whose formal accuracy is higher than two (away from the isolated points) using 
a minimal band width, was presented in [21]. 

In the present paper we are considering two—level, (and for simplicity 
only), explicit finite difference approximations to a scalar conservative law 
having two properties: 

(1) Lax-Wendroff conservative form [15], and 

(2) Roe's incremental form [22], 

In Section 1, we introduce a modified flux and a viscosity coefficient, 
both grid dependent quantities, with which we set up and advocate a third form 
of these schemes the viscosity form. A sufficient TVD criterion is then 
presented in terms of this form. 
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In Section 2 we consider schemes which have a 3-point stencil. In this 
case the modified flux coincides with the original (differential equation's) 
flux; accuracy is consequently limited to first-order. 

Thus, in Section 3, we consider wider stencils. A comprehensive framework 
of second-order resolution (SOR), total variation diminishing (TVD) schemes is 
presented. The viscous form of the scheme again plays a key role here. 

Two concrete examples for the construction of SOR-TVD schemes according to 
the above guidelines are detailed in section 4; both approaches use the 
viscous form of the underlying 3 point TVD schemes. Special attention is 
called to the second ("piecewise constant viscosity modification") recipe in 
Corollary (4.9), which results in an easily implemented and highly attractive 
two step formulation. 

These four sections comprised Part I. In our attempt at a unified 
treatment of this subject, we derived (and in some cases rederived) some 
notable conclusions: 

(i) Three point TVD schemes are, at most, first-order accurate. 

(ii) Accuracy at nonsonic critical points is limited to first-order, 

(iii) Our piecewise linear flux correction (Theorem 4.3) extends Harten's 

recipe [10], under a CFL of limitation 2/3, (Remark 4.9), in 
agreement with that found by Sweby in [24], 

(iv) Our recipe for converting three point TVD schemes to five point SOR- 
TVD schemes seems quite general and attractive from both a 
computational and a theoretical point of view. The latter point of 
view will be used in Part II. 

The discussion concerning a cell entropy inequality begins in Section 5 
using the canonical Godunov scheme. Section 6 contains the heart of our 



discussion. The grid dependent modified flux of Section 1 is extended to a 
modified flux function defined on the intervals connecting the grid values . 
The definition involves a piecewise linear function (i.e., the "double wing" 
seen in figures (6.4)) and the remark made after the statement of Theorem 
(6.11), which takes into account the presence or absence of critical points. 

The schemes under consideration are then expressed as convex combinations 
of Riemann problem solvers, i.e., of Godunov type schemes, using the above 
modified flux functions. We obtain an entropy in cell inequality involving a 
familiar (integral) residual term. 

Section 7 discusses the significance of E fluxes (introduced in [18], 
c.f also [26]). We show that for any approximate entropy in cell inequality 
to be valid for all intervals, the underlying scheme must have an E flux 
hence be at most first-order accurate (Theorem (7.2)). Conversely, the 
existence of an E flux implies the nonpositivity of the residual term 
mentioned above, hence implies a general cell entropy inequality. To obtain 
convergence for SOR schemes we give up the requirement that all discrete 
entropy inequalities be valid. This limits our convergence proof to convex 
(or concave) conservation laws, as in [19], [20], 

In Section 8 we treat SOR-TVD schemes satisfying a single quadratic cell 
entropy inequality. Here we use the viscosity modification recipe of Section 
4. A cell entropy inequality follows by estimating the residual integral 
presented in Section 6. The rather delicate tuning of the modified flux so as 
to comply with the three criteria SOR, TVD, and entropy inequality, leads us 
to a wide class of explicit in time SOR-TVD schemes (again only for the convex 
or concave scalar case). 

The main convergence results are contained in Theorems 8.4 and 8.5 
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PART I. TOTAL VARIATION DIMINISHING SCHEMES 


1. The Modified Flux and the Numerical Viscosity Coefficient 
We study two-step difference schemes of the form 


v v (t+k) = H(v v _ p (t),...,v v+p (t); f,x), 0.1) 

serving as consistent approximations to the scalar conservation law 


3u 

3t 


&-<*■ 


( 1 . 2 ) 


Here, v^(t) = v(x^,t) denotes the approximation value at the grid point 
(x = vAx,t), k = At, and Ax are respectively, the temporal and spatial mesh 
size with fixed mesh ratio X = k/Ax, and p a natural number. 

We postulate 


Assumption 1.1 

(i) The scheme (1.1) admits a conservative form in the sense of Lax- 
Wendroff [15] 

H (Vp- — 'Vp ! f,l) ■ \ ■ x C h v + V 2 ‘ Vv 2 ) <U3 > 


where h^ + y^ stands for the Lipschitz continuous numerical flux 


h i. = hfv .!,••• »v . ; f ,x), 
v +72 v v-p+1 v+p '* 


consistent with the differential one 
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h(w,w,»»»,w; f,X) = f(w). 

(ii) The scheme (1.1) can be also written in an incremental form 

Hfv ,»»»,v ; f,X) = v + C + i, Av 1 - - C i. Av i. ; (1.4) 

v v-p v+p 1 v v+ V2 v+ V2 v- V2 v- 72 

where we use the standard notation 


Av , i, = v , , - v . 
v+ x /2 v+1 v 

Equating the right hand sides of (1.3) and (1.4), rearranging and dividing 
by X, we find 


V v 2 + r C i / 2 4 v v 2 " V V t V v 2 4 v v 2 • <1,5a) 

We term the equated grid dependent quantities in (1.5a) as the modified flux 
associated with scheme (1.1): 


g v h v± V2 + X C v± V2 A V V2 

■ i t h v- v 2 + V1/2 + r V v 2 4 v v 2 + 1C v 2 v 2 J • 

As in [25, Section 2], we now use the consistency relation 

c- v - c + ... x ly 2 , 
v+V 2 V+V 2 Av v +1/ 2 


(1.6a) 
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to concludee that with a given modified flux, g^, there is only one degree of 
freedom in setting up the recipe of the difference scheme (1.1): this will be 
manifested in terms of the quantity. 


^v+ V2 C v+ V2 + C v+ V2 * 

Indeed, by averaging (1.3) and (1.4) we find 


(1.6b) 


v v (t+k) - v v (t) - | (h v+ l/ 2 - W 2 ) + 7 (Cv 2 A W 2 ‘ 

adding and subtracting ^ ( c v +l/ 2 Av v +l/ 2 ” C y-l/ 2 Av v -l/ 2 ) on the right-hand 
side, we arrive at 

V v (tHO = v v (t) [(Vv 2 + lVi/2 i W2 ) " ( \- 1 / 2 + T < V 1 4 4 VV 2 )] 

+ 7 [ cV V 2 + < + v 2 ) 4 V V 2 - (V V 2 + <- V 2 ) 4 V V 2 ] • 

In view of (1.5b) and (1.6b), the scheme (1.1) is thus finally recast into its 
viscous form 


V t+k) " v v (t) "T (g v +i ■ Vl) + T (1 * 7a) 


expressed in terras of the numerical viscosity coefficient [10,26] 

g y * g y+l ~ 2h v+V? 


Q v+ i /2 = c v+ i / 2+ c v+i/2 = X- 


W 2 


(1.7b) 



Let TV[v(t) ] = £ v |v v+1 (t) ~ v v (t)| denote the total variation of the 
computed solution at time level t; the following lemma provides us with 
sufficient criterion for the scheme (1.1) to be total-variation diminishing 
(TVD), in the spirit of [10]. 


LEMMA 1.2 ( Total-variation diminishing) 

The scheme (1.1) has a diminishing total-variation 


TV[v(t+k)] ^ TV[v(t)], 


provided its numerical viscosity coefficient Q v+ \j , satisfies 

.K+isJ 


l VV 2 


= ^V+ V2 = 


( 1 . 8 ) 


Proof . By averaging (1.6a) and (1.6b) we find 


cU /2 


In view of these last relations, the inequalities (1.8) boil down to 


>°> C „ + lA^°. 


v+ V 2 = °* W 2 ^ 


and TVD follows along the lines of [10] 
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2. Three-Point TVD Schemes 

In the case of 3-point schemes, p = 1, it was shown in [25; Lemma 2.1] 
that there exists one and only one incremental form (1.4), whose incremental 
coefficients are 


f - h . I, 

(T +1/ ■x-g - . -S t V 2-, C-J.-X 

v+V2 A VV2 V /2 


f v+1 *\>+ V? 

A vv 2 


( 2 . 1 ) 


Inserted into (1.5b) we find that the modified flux in this case coincides 


v v (t+k) = v v (t) (f v+1 - f v-1 ) +| [A(Q v _l /2 Av v _ 1/2 )]. (2.2) 


Abbreviating a(w) = -y— (w), the following notation will be used throughout 


W 2 



‘vv 0 

(2.3) 

A VV 2 = °* 


The TVD constraint (1.8) now reads 

x K+i / 2 I ^ w 2 - 


(2.4) 


and the following list quotes the most frequently referred to schemes in this 



-10- 


3-point TVD category [6,7,16,22]: 


Cv 2 - 


(Lax-Friedrichs scheme) (2.5a) 


0^+1/ = . * — / |a(v)|dv, (Engquist-Osher scheme) (2.5b) 

2 v+ V 2 v 


a G _ Max 

VV2 X (v-v y ).(v-v v+ i) < 0 ' 


l W 2 


(Godunov scheme) (2.5c) 


&V 2 - x K+V 2 I- 


(Roe-Murman scheme) (2.5d) 


In fact, in the special case under consideration of 3-point schemes, condition 
(2.4) is necessary as well as sufficient for TVD, see [25, Corollary 2.3]. 
Hence, the following 3-point schemes are not TVD ones 


q£ v = (* a v+ l/ 2 ) 2 > (Lax-Wendroff scheme) (2.6a) 


Cv 2 -°- 


(Forward-Euler scheme). (2.6b) 


A special significance is attached to the Lax-Wendroff scheme 


v v (t+k) = v v (t) - | (f v+1 - f v+1 ) + ± [A(Ua v _ 1/z ) 2 Av v _i /2 )]; (2.7) 


when seeking second-order accurate schemes, the further limitation placed on 
Q^ + in the 3-point case, singles out the Lax-Wendroff choice, (2.6a). In 
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view of the above, we therefore conclude that 3-point TVD schemes are at most 
first-order accurate [10, 18, 25]. 

At this point we are widening our discussion to include schemes whose 
stencil occupies more than three points: the further freedom in setting up 
the modified flux g in this case, would enable us to achieve higher (than 
one) degree of accuracy. 


3. Second-Order Resolution Schemes 

We start by rewriting the modified flux, g^, in terms of the correction 
g^ to the original one, 


1 - 


(3.1) 


Inserted into (1.7a), our scheme now reads 


v v (t +k) = v v (t) - 2 [(f v+1 + x I v+1 ) " (f v _! + x 

+ 7 [ a(q v-V 2 a VV 2 ^ 

= v v (t) " 7 ( f v+l ' f v-l) + 7 [ A <V V 2 A V V 2 )]» (3 ‘ 2) 

where Q^ + stands for the modified viscosity coefficient given by 

*V + Vl 1 m J f v + f v+l 2h v+ V 2 

a vv 2 L a vv 2 


W Vv 2 


(3.3) 
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In other words, we rewrite our scheme w.r.t the original flux , f^, modifying 
the viscosity coefficient instead. 


Remark 3,1. The viscosity coefficient discussion in [10], corresponds to the 
modified one in (3.3); indeed, in the 3-point case treated in [25] where 
g^ = X(g^ “ f y ) = 0, it coincides with (regular) viscosity in (1.7b). 

Comparison of (3.2) and the Lax-Wendroff scheme (2.7) leads us to the 
following second-order accuracy requirement 

S., "b S.j-1.1 9 9 

%+ V 2 5 V V 2 av v+ i / 2 = x (a \+ V 2 ; (3,4a) 

the deliberately vague notation of X^(a^)^ + i^ on the right, stands for 

Q LW i. = X 2 (a , i, ) 2 modulo first-order errors, i.e., 
v+ V 2 v+ V 2 

i 2 (> 2 ) v+ l/ 2 = ^ 2 (v l/ 2 ) 2 + °< M),, y 2 . C3.4b) 

Together with the TVD constraint (1.8), we finally arrive at a general 
description of second-order TVD schemes, which is summarized in the following 

LEMMA 3.2 (Second-order TVD schemes) . 

Consider the difference scheme (1.1) written in its viscous form 


v v (c+k) - v v (t) - | (g v+1 - g^j) + i [A(V 1 , 2 


with modified 


flux , g^, and viscosity coefficient , Q^ + 


by 


Av v _l /2 )], (3.5a) 

, given respectively 
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22 s v+l 

Q v+ i/ 2 - * < a \u. 1/. + * 


'\>+ V 2 Av 


v+ V 2 


(3.5b) 

(3.5c) 


The scheme (3.5) Is second-order TVD, provided the correction terms g and 
2 2 

X (a )^ + on the right-hand side of (3.5b-c) satisfy the following two 

requirements of 

(i) second-order accuracy: 


* 2 <» 2 \, + y 2 - (»v y 2 f + «( |av |, /2 ; 


(3.6a) 


(ii) total-variation diminishing 


” +1/ 2 4v „+1/ 2 


2 2 

S X (a ) v+l/ 2 + ” 


V1/2 


< 1 . 


(3.6b) 


As we shall see below, one cannot satisfy both requirements (3.6a) and (3.6b) 
at the nonsonic critical grid values v^ where Av^_ y • Av^ + < 0 * a(v^), 
and, therefore, second-order accuracy must be given up at these values. 
Difference schemes with (formal) second-order accuracy at all but those 
critical grid values are classified as having second-order resolution, after 
Harten [10]: thanks to the TVD property, no new such first-order accurate 

critical grid values are added during the computation (since the scheme is 
monotonicity preserving [10,25]), and the overall second-order accuracy does 
not seem to be degraded in this case, at least in the norm. 

There are various approaches to the construction of second-order 
resolution (abbreviated hereafter SOR) schemes, the main three of which are 
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the following: 

(i) The MUSCL approach [27,19], directly modifies the first-order 
numerical flux h^ + y^ , using a Godunov-like solver for a limited slope 
piecewise linear grid data. 

(ii) The modified equation approach [10,11] directly modifies the 
flux fCv^), based on considerations of the modified equation associated with 
a first-order scheme, so that the limited modified flux constructed, 
guarantees second-order TVD resolution. 

(iii) The flux limiter approach [24] directly modifies the first-order 
incremental coefficients C* + y^ , using a class of flux limiters which 
preserve TVD and maintain second-order resolution; in particular, such flux 
limiters can be chosen to interpret MUSCL and modified flux-type schemes due 
to Van Leer, Roe, Harten, Chakravarthy, Osher and others, see [24], 

All the above approaches can be entertained of course, within the general 
framework provided in Lemma 3.2. A common "limiting" feature shared by all of 
these approaches can be directly derived from the TVD constraints in that 
lemma: the grid values g y read in (3.5), are expected to form a first-order 
"smooth" grid function correction, at least at the generic noncritical zones; 

2 

since (a ) y+ 1/^ is determined up to first-order perturbations, see (3.6a), 
we can absorb such a perturbation, ±Ag y+ y^ /Av^ y^ = 0( |Av| y^ ), into the 
middle term of (3.6b), and the right-hand inequality now reads 


, i., , o v+1 

x u W/ 2 + 2 


VVa +2 Sv^‘- (3.7a) 


The left-hand inequality in (3.6b) gives us 
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3 v+ V2 


' X (a") 


V 2 * 2 i 


V v+ V 2 


X VV 2 " X 2 (a 2 ). 


v+ V 2 


vv 2 


(3.7b) 


We shift indices in the first and fourth inequalities in (3.7). The resulting 
inequalities boil down to the following inequality to be satisfied by TVD 
schemes: 


1 («v l/ 2 - ‘V’vi 1/2 ) 5 - i J (1 - X V 2 )• 

Hence, g^/Av^jy^ must lie between the above two bounds, both determined up 
to first-order perturbations . In particular, in view of the second-order 
accuracy requirement in (3.6a), the sign of the parenthesis on the left of 
(3.7c) is determined by that of the first terra, ±Xa y± ; if further, we are 
at a nonsonic value, a(v^) * 0 , then this first term cannot be absorbed as a 
first-order perturbation of the second one, and (3.7c) requires the middle 
terms 8 v / Av v ±l /2 to be positive. This implies that the second-order accuracy 
must be given up at the nonsonic critical values v where 
Av v _ 1 / 2 » Av v+ iy, 2 < 0 * a(v v ), [25,20]. 

We now describe two specific recipes which convert arbitrary 3-point TVD 
schemes into second-order resolution ones, in the spirit of the above 
guidelines. Indeed, although the reasoning may be different, the various 
"conversion" recipes including the two below, end up with difference schemes 
which bear close similarities to each other, as dictated by the framework 
provided in Lemma 3.2. 
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4. Two Recipes for the Construction of SOR-TVD Schemes 

We start by Identifying a 3-point TVD scheme by its numerical viscosity 

coefficient Q , i, = Q(f:v ,v , ,X): it is considered here as a functional in 

N+ V 2 v ’ v+1 ’ - 

the original flux f(») with the further dependence on the grid values 
v v > v v+ l anc * t * ie mes h ratio X. In view of the necessity of the TVD 
constraint, (2.4), we may assume without loss of generality that the following 
functional inequality is obeyed: 

Iz(w-) - z(w )| 

X | » 2 - “1 I - qU; W X) ' (4a) 

Thus, the TVD requirement (2.4) is reduced in this case to the CFL-like 
condition read in the right of (2.4) 


vv 2 


(f) = Q(f;v v 1) 5 1. 


Remark 4.1. The above description seems to exclude several difference 

schemes, where it is actually the inequality on the left of (2.4), 

XIAf , 1 , /Av 1 . I < Q , 1 , (f), which leads to the CFL limitation, as in the 

1 v+ 72 v+ 72 = v+ 72 

case, for example, of Lax-Friedrichs scheme (2.5a) with, ^ = 1. Together 
with the TVD constraint (2.4), however, the Lax-Friedrichs scheme, for 

example, can be equivalently represented by a numerical viscosity coefficient 
Q^ + = Hax(l,x|Af v+ ^ |/|Av v+1 ^ |) which fits into our above interpretation. 

In general, an expansion of Q(»;*,*,X) in powers of X is called for, but 
we shall not elaborate on that here. 

To construct an SOR-TVD scheme, the two first-order correction terms g 

^v+ V2 
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relation (3.5c), we shall also make use of the further (first-order) 
flexibility in determining the viscosity coefficient, • 

Let us denote 


■v = iK-v 2 + vi/2 5, w 2 5 sgn ( 4 vv 2 L <4 - 3a) 

and, in view of (3.7c), set the flux correction to be 

Min ([ < *v±v 2 <£) - x2(a2) v±v 2 i • K±v 2 d- < 4 - 3b > 


2 2 

Here, the term X (a ) i, is chosen so that 
v+ Vo 


X 2 (a 2 ) v+V 2 -(Xa v+1/2 ) 2 + 0(|Av|) v+1/2 , »V )<t ^ S V ^ ( £ ) (4.3=) 


with otherwise arbitrary first-order perturbation 0(|Av|)1^ • The modified 
flux is then given by 


(4.3d) 


and we are considering the difference scheme whose viscous form is expressed 
in terms of that modified flux 


v v (t+k) = v v (t) (g v+1 - g v-1 ) + \ [A(Q^ R i / 2 Av v _ 1/2 )]. (4.3e) 


We have 
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THEOREM 4.2. 

The difference scheme (4.3) is SOR-TVD provided Its viscosity coefficient , 

CV 2 • satlsfles 


a , i -aJfr 

v+1/ 2 *VV 2 


VV 2 


l Q v+l/ 2 - W 2 (f)| * 0<|Sv|) v+V 


(4.4b) 


Remark . We note that the definition of (a in (4.3c) still allows the 

flexibility of (some limited amount of) first-order perturbation as long as 
2 2 

Q^ + y^ (f) - X (a )^ + y^ remains nonnegative; the latter guarantees that the 
sign of g y agrees with s y± y^ . 


Proof . The first condition, (4.4a) is nothing but the TVD requirement (3.6b); 
consult (1.8). Away from the critical values where s^ = 0, (4.4b) implies 
second-order accuracy; consult (3.4), 


SOR _ s y + g y+l 
v+1/ 2 Av v + 1/ 2 


- <C\- fvv 2 (f) - x2(a \> + V 2 l + ° 0 av|) v+ i /2 

= [Q^+v 2 - W 2 (f) ] + x 2( W 2 ^ 2 + 0( l Av l ) v+1/ 2 


W 2 J 


V v 2 • 


Two concrete choices for SOR viscosity will now be discussed. 

As a building block for the first, we introduce the piecewise linear flux 
correction 
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**+v 2 <T > V ♦», 

and use the original 3-point viscosity functional, setting 


(4.5a) 


C / 2 ■ Q t f(,) V Vi- »)• (4 - 5b) 


Lipschitz continuity of the viscosity functional implies 


( 1 ) 


ICvwH s . 


and we end up with 

THEOREM 4.3. (Piecewise linear flux modification) 

The difference scheme (4.3), (4.5) is SOR-TVD under the CFL-like 

condition . 

VV 2 (£) i ITT' (4 - 7) 


Proof . The SOR requirement in (4.4b) is fulfilled because of (4.6), whose 
right-hand side is of order 0(|Av|) v+ ^; the TVD requirement on the left of 
(4.4a) follows from (4.1), and we are left with the inequality on the right, 
requiring 


To be precise, the constant L stands for the Lipschitz constant times 
the (assumed finite) maximal bound of the ratios 

IKW 0 ' v T v+r 
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^v+ V 2 


(4.8) 


Indeed, in view of (4.6) we have 


Cv 2 * Vv 2 <£> + L 


A vv, 


w 2 


(4.9a) 


since and § v cannot have different signs by the choice of s in 

(4.3a), we also find 


Ag v+ l h 


W 2 


Max (|g I, |g J) 1 2 2 

* - | - 4V +1/ J S 2 [V Vj (£ > - 1 <* V V 2 ]■ (4 ' 9b > 


The last two inequalities yield the crude bound 


CV 2 < V l / 2 (£) + \ (V l / 2 (f) - »V >*. V 2 ] S ( 1 + |)V v 2 <£) < 4 - 10) 

and augmented with the assumed CFL-like condition (4.7), the TVD constraint 
(4.8) is now satisfied. 


Remark 4.4 . With the primary examples of Engquist-Osher, Godunov and Roe- 
Murman schemes in (2.5b - d), the Lipschitz constant L in (4.6) may be taken 
to be L = 1, and the CFL limit in (4.7) is found to equal (compare, e.g., 
[24, Section 3]) 

Vv z (£) ^ 2/3 - <4 - u) 


2 2 

(Taking into account the negative term -X (a ) 


v+ V 2 


which was ignored in the 


crude bound of (4.10), the CFL limit may be further improved to be 
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W f) 


2 + *V> v+ l/ ? 


(4.12) 


Remark 4.5. In the case of Lax-Friedrichs scheme, where 0^l, - 1 violates 

- v+ 72 

(4.11), we may take L to be zero, leaving us with the original 3-point CFL 
limitation 


»|a v+ l /2 l $ I- (4.13) 

Remark 4.6 . Consider the special choice (a^)^ + j j^= (ain (4.3b), 
which is admissible in view of (2.4). The piecewise linear flux modification 
in Theorem 4.3, extends Harten's recipe [10], which was restricted to 
viscosity functions of the form Q^ + y = |); (hence, the Engquist- 
Osher and Godunov schemes, for example, were excluded from the discussion 
[ 10 ]). 

We now turn to a second more attractive choice of SOR viscosity, inspired 
by Theorem 4.2 (see also [10, p. 368]). Here, we use a piecewise constant 
viscosity modification of the form 

Cv 2 - W f) * 


v v? I 


W 2 | 


(4.14) 


which leads to 

THEOREM 4.7. ( Piecewise constant viscosity modification) 

The difference scheme (4.3), (4.14) is SOR-TVD under the CFL-like 

condition 

Vv 2 (f) s 2/3 - 


(4.15) 
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Proof. The SOR requirement in (4.4b) is directly verified by the very 
definition of y in (4.14); the TVD requirement on the left of (4.4a) 
follows from (4.1), 



< X 

f(v v+l> - f <V 


4 VV ? 

Xa , 

v+V 2 Av v+1 /“1 

V v+1 " \ 


4v v+1/ 2 


✓ n S0R 

s W 


(4.16) 


and we are left with the inequality on the right, requiring 


Vv 2 


(4.17) 


Indeed, the same estimate used in Theorem 4.3 before, see (4.9b), together 
with the definition of , (4.14) gives us the pessimistic bound: 

Cv 2 S V 1/2 <£> + i ['V v 2 (£) - 1 / 2 ] i 1 V i / 2 (£) (4a8) 


and (4.17) follows in view of the assumed CFL—like condition (4.15). 


Remark 4.8 . As before, by improving the crude bound of (4.18), the CFL-limit 
may be improved to be that of (4.12). 


We note that the piecewise constant viscosity modification involves linear 
corrections of the numerical flux f in (4.3d) and the viscosity coefficient 
in (4.14); both are appearing linearly in the original 3-point scheme we are 
starting with; see (2.2). Hence, the resulting modification recipe does not 
change the underlying 3-point TVD code, but, rather, adds to it an 


antidiffusive term which boils down to 
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-1 OUi - Vi) + 4 v v 2 1 


(4.19) 


To make our point more precise, we state as our final result of this first 
part, the following 

COROLLARY 4.9. ( Piecewise constant viscosity modification - revisited) . 
Consider the 3-point scheme 


v v (t+k) - H( Vl , v Vl ; f,x) S v v - x(h v+1/2 - h v _ 1/2 ) 


(4.20a) 


and assume the TVD-like constraint, compare (2.4) 


x i v i /2 1 i v v 2 (f> i I • v v 2 -» fv 1 • (4 - 2ob) 


holds. Then the modified scheme 


v u (t+k) - H(v v _ r v Vl ; f,X) - 




(4.21a) 


is SOR-TVD. Here , g y is the flux correction, given by 

2 , 2 . 


i y -^ M in ([Q^^Cf) - ^^ )v±1/2 ] |Av v±1/2 |) ( 4 . 21 b) 


where 


* < a W 2 =X 2 K +1/2 ) 2 + 0 (|Av |) v+1/2 < Q v + l / 2 (£)• 


(4.21c) 



- 24 - 


Alternatively, we can rewrite (4.21a) in the conservative form 


v v (t+k) = H(v v _ 1 , v v> v v+l ; f,x) - [K v+ i / 2 - E v . 1/2 ] (4.22a) 

- - x[(h v+ 1 / 2 + |K v+1/2 ) - (VV 2 + T K v-V 2 )I - 


SHi Vv 2 


is the numerical flux correction given by 


v+ V2 


■is 


V+ V2 ^v+1 ” I 


(4.22b) 


Remark 4.10. The Lax-Friedrichs scheme violates the CFL limitation (4.20b). 
Nonetheless, an "antidiffusive" modification follows in this case directly 
from the flux modification of Theorem 4.3, with numerical flux correction 

\>+ V 2 = 2 ( + ®v+l^ » 


Remark 4.11 . The conservative form of scheme (4.22) is identical with the 
corrective type of Harten's artificial compression method [8,9], the 
difference lying, of course, in the exact details of the flux correction used. 
In particular, we have an operator splitting which could be easily implemented 
as a two-step predictor-corrector method [8, Section 6]. 
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PART II. ON A CELL ENTROPY INEQUALITY 

5. Godunov's Scheme and its Cell Entropy Inequality 

An entropy pair (U,F) associated with the conservative model (1.2), 
consists of convex entropy function , U(»), augmented with an entropy flux , 
F(•), such that U'f' = F'. The requirement of having an entropy inequality of 
the form 

(u(x,t)) +|| (u(x,t)) * 0 (5.1) 

for all entropy pairs, singles out the unique physically relevant (weak) 
solution of (1.2), e.g., [13,14]. Accordingly, we are seeking conditions 

which guarantee that the difference scheme (1.1) will satisfy a cell entropy 
inequality of the form [12] 


U(v v (t+k» < U(v v <t)) - *(F v+ 1 /2 - VV 2 )S <5-2s) 

here, F v+ 1/^ = F ( v v _q+i»* * *» v v +q» *- s a numerical entropy flux, assumed 

to be consistent with the differential one^^ 

F(w,w,«*«,w; f,X) = F(w). (5.2b) 

If this holds for wide enough class of entropy functions, the cell entropy 
inequalities (5.2) will guarantee the uniqueness of limits of all converging 


^ ^ Both the differential and the numerical entropy fluxes are denoted by F; 
their distinction is made by the number of arguments. 
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subsequences of the total-variation bounded solutions, and there follows the 
convergence of the whole computed sequence to the unique physically relevant 
(weak) solution , e.g., [3, 12, 18, 23, 26]. 

As an example, we begin with the all important cell entropy inequality 
associated with Godunovas scheme : by averaging two Riemann solvers, one 
obtains Godunov's scheme whose numerical flux is given by [17,18] 


h G i, = s . i, Min fs i. f(v)l, s i. = sgn(Av . i, ). (5.3) 

v+V 2 v+V 2 ( ).( v -v ) < 0 v+1/z V /2 ^ k 

V v +1 

Assuming that the minimum on the right-hand side is attained at a (not 
necessarily unique) value, v G + = v G (v v ,v y+1 ,f), 

s »+Vi, < o [ VVi ,<v) ] ■ £( vV 2 >■ <5 - 4a) 


then the Godunov scheme 


v(t+k) = H G (v v . 1( v v ,v v+1 ; f.X) 


(5.4b) 


boils down to 

hG (Vi* v v>Vi ; f » x ) = v v " x [ f (VV 2 ) " f Cvv 2 ^* <5,4c) 

Integrating over a typical cell the differential entropy inequality (5.1), 
then Jensen's inequality yields the following well-known lemma, e.g., 
[3,18,26] . 
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LEMMA 5.1. ( Cell entropy inequality of Godunov scheme .) 

Consider the Godunov scheme (5.4) satisfying the CFL condition. 

X Max |f'(v)| < 1. (1) (5.5a) 

Then, the following cell entropy inequality holds 

u[h G (Vi’VVi ; f ’ x )l ^ u(v v ) ‘ x [ f( V 1 / 2 ) " f( VV 2 )] (5,5b) 

for all entropy pairs (U,F). Equivalently, (5,5b) can be rewritten 

G 

G VV 2 

U [ H ( v v _l’VVl ; f,X ^ = U(v v ) " / U'(w)f'(w)dw. (5.5c) 

VV 2 

As in [26], the cell entropy inequality associated with the Godunov scheme, 
will be used as a building block for studying such an inequality in 
conjunction with other TVD schemes. To this end, we shall make use of a 
modified flux function introduced below. 


The maximum is taken over all values v varying between v v _^v v , and 

v V+l* 
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6. The Modified Flux Function 

While studying total-variation diminishing schemes, only the modified flux 
grid values, g y , played a role, see the TVD requirement (1.8); subsequently, 
a simple piecewise linear modified flux correction, g v+ iy^(v), was introduced 
in (4.5a), in connection with the (first) recipe of SOR-TVD schemes. 

Regarding the question of a cell entropy inequality which we now consider, 
a more sophisticated construction of a modified flux function , is required: 
{r (v+ / 2 )( v j denote its restriction to the I y+ “ interval. 


let g v 


T v+ l/ 2 “ { v l ( v-v v )( v-v v+ i) £ 0}. T he piecewise linear function g^ V+ ^ ) 


(v) 


interval, through a constant numerical flux value, h v+ l/^t it depends on two, 
yet to be determined parameters (s + ,s ), s* > 0, in the following fashion: 
Setting the intermediate values 


vv 2 
vv 2 

we define, depending on whether v y ^ v y+ j 
fg v - s + (v-v y ) 


+ , + x . 1 
2 W 2 (s ) -v ^-7 

5 VV 2 (s " ) ‘ Vi •; 


(V+V 2 ), + -x 

g 7 (v; s ,s ) 


h v+ V 2 


K ' Vv 2 ) 

r ( g v+l “ h v+ V 2 ^’ 

or V V+1 ^ V 

T v = v^ V 1/2 - 

+ - 
v v+ V 2 - V = V V+ V 2 

v; + i /2 < v < v v+1 


(6.1a) 

(6.1b) 




(6.1c) 
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Vl * V * VV 2 


g (v+1/2) (v ; s\s‘> - hv+1/2 , V 1/2 , v v+1 < V,. 


’ W 2 S V SV 


Assume the TVD condition, (1.8), holds 


‘ Av . i, v+1 


Then g^ V+ ^ \v; s + ,s ) given in (6.1) and satisfying 


is a well-defined piecewise linear function in the I i, interval. 


Verification. In view of (6.2a) we have, using (1.10), 


C \>+ V 2 = °* C v+V2 = ^v+ V 2 ” C v+V2 + C v+V2 - ** 


and the definition of g in (1.5b) yields 


V1/2 ■ v v + T C g v - h v+ V 2 5 ■ % + 7 + C v+ v 2 4 v v 2 
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Figure 6.4a 





where the graph of 


A completely analogous situation occurs when V v+ j < v v » 
g(v+ I 2 takes the form, see Figure (6.4b). 



Our next lemma is in the heart of the matter. 


LEMMA 6.2. 

Assume the TVD condition, compare (1.8 ), 


S <Uu. 


< 0 < 1 
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HG [VVVr g(v+1/2> ( s+ ‘ s " " ex_1 )> 0 ' X 1 " v v ' 0 ‘ lx tVv 2 ' s vl 


^V,.VV* (rl41 ( st ■ s ' ■ ex " 1 )- 9 ' lx l ■ \ - •“‘I*. - V v 2 1 


Proof . We first note that according to Lemma 6.1, g^ V± ^ )( s + = s = ex *) 
are well defined: (6.5) yields 


By definition we have, see (5.4c), 


k g [v ,V V ;g (v+1/ 2V - s- = ex' 1 ), e'h] 


■ \ ~ e" 1 Xtg <v+l/2) (»^ ( .i /2 ) - g ( ' l ' 1/2) (v v )]; (6.7a) 


straightforward computation gives us, see (5.4a), 


g (v+ V 2 )(G ) = s y Min [s y g (v+ V 2 >(v)] = h y (6.7b) 


W 2 J = W 2 ““ L W 2 g 

V+ V 2 


while consistency implies 


g (v+1/2, (v ) 


and (6.6a) follows. Similar arguments apply for (6.6b). 
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Remark 6.3. Godunov's numerical flux is determined by the differential one, 
computed at the intermediate value v ^\j ; see (5.4a). In our case, (6.7b), 
any intermediate value between v^ and v^ ^ can be taken: for later 
purposes, we shall choose in particular 


V>4 (,) "?tvv 2 (s+ 


W 


• ex' 1 )] 


_s±L —ir, -,i. 

20 ^u+l 8 u7* 


( 6 . 8 ) 


Taking 0 = 1/2 in Lemma 6.2 and averaging (6.6a) and (6.6b), we finally 
arrive at 


THEOREM 6.4. 

Consider the difference scheme (1.1) given in its viscous form (1.7a) and 
satisfying the TVD condition 

i V ‘/ 2 * 1/2 - (6 - 9) 


K* 1 /? 

i 4 w 2 


We then have 


\< t+k > - V - s ' ■ it). 2 M 

+ hG [ V v> V v>Vi ;S<V+ 1/2> ( S+ ‘ s " ’ 2> ]]- (6-10) 

In other words, we have shown that any TVD satisfying difference scheme, 
(1.7a), (6.9), is given as an average of two Godunov solvers, thus refining a 



-34- 


similar (convex) decomposition introduced in [26, Theorem 5.1]; the current 
decomposition differs however, in its use of a modified flux function, rather 
than using the original differential one. In fact, away from the critical 
extremum values v y where Av y _ • Av^ + <0, we can do even better: 
under appropriate TVD condition, any difference scheme coincides with one 
(rather than the average of two) Godunov solver. This is the content of our 
next 


THEOREM 6.5. 

Consider the difference scheme (1.1) given in its viscous form (1.7a) and 
satisfying the TVD condition 


X 


Ag **V2 

Av v±V 2 


< Q, 


! v± V 2 


( 6 . 11 ) 


Then, in the neighborhood of noncritical values v such that 
Av v _V 2 * Av v _ 1/z > 0 , w e h ave 


v y (t+k) = H G [v^_^» v v , v v+1 5 g(s + = s" = j), X]. (6.12a) 


Remark 6.6 . The modified flux function g(v) = g(v; s + ,s~) appearing in 
(6.12a), is uniquely defined in this case, composed of its restriction to the 
consecutive intervals I^ ± , see Figures (6.13), 


g(v; s + = s“ = -J-) 


J ( V" V2 )^ v . , + „ . j). 

V £ V >/ 2 

, ( ' ,+ 1/ 2 > (v, S + - -{), 

v s VV 2 


(6.12b) 




Figure (6.13a) 


g(v; s" = s" = i) 



Verification of Theorem 6.5 is straightforward: by (6.11) we may employ Lemma 
6.2 with 0 = 1 in I v± ly^ » using (5.4c) and (6.7b) the result (6.12) 


follows. 
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We next use the cell entropy inequality associated with Godunov's scheme, 
see Lemma 5.1, to infer such an inequality for the difference schemes under 
consideration. We start with 


LEMMA 6.7. 

Consider the TVD scheme (1.7a), (6.9). Then for all entropy pairs 

(U, F) we have 


U[v v (t+k)] i U[v v (t)J - F v _ 1/2 ) 

V 

+ A / U"(v)[g^ V ^(v; s + = s = 77 ) “ f(v)]dv 

G 

VV 2 

VV 

+ A £ U"(v)[g^ V+ ^(v; s + = s = “ f(v)]dv, (6.14a) 

with a consistent numerical entropy flux given by 


F v+ i/2 = f( VV 2 ) " U '(VV 2 ] * [VV 2 ■ f K + v 2 )]- (6 - 14b) 


Remark . Here, the intermediate Godunov value v^ + may be taken as, see 

( 6 . 8 ), 

V l /2 = vj + i /2 (9 = V 2 ) = - ■ 2 - 1 - “ J (§ v+1 " g v )« (6.14c) 

Proof . The maximal slope of g^ V+ ^ ^(v; s + = s = 7^) is V 2 A, 
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0 <v±V 2 r (v . / ^ 


Therefore, we may apply Lemma 5.1, obtaining, see (4.5c), 


u [ h< H-i’VV s ( v_1/2) ( v 5 s+ = s " - 2T^» 


< U(v) -2X / U'(v)g^ v ^ (v; s + = s = -orOdv (6.15a) 


; g< V+1 / 2 )( v; s + = 


^ U(v ) - 2X / U'(v)g^ V+ ^ 2 ^ (v; s + = s = -^r-Jdv. (6.15b) 


Adding and subtracting 2X[F(v^) - F(v^_ )] =2X / U' f' to the right- 

vv 2 

hand side of (6.15a), we find after integration by parts 


d [“ G [Vi*VV B+ ' s " ■ it 2 >l] 


i U(v v ) - 2X t* 0'(g (v - V 2 > - f)' - 2X[F(V V ) - F(v^. 1/2 )] 

VV 2 

v=v 

= U(v v ) - 2X • U'(v) • (g (v ~ 1/2> (v) - f(v)) V - 2x[F(v v ) - 

V= vv 2 

v 

+ 2X / U"[g (v-1/2) - f]. 

G 

VV 2 


(6.16a) 
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In a similar manner, (6.15b) yields 


u[h G [v 


(V+V2 ) 


(v; 


■ 2X 


r )• 2 x l] 


< U(v v ) - 2X . U'(v) . (g <u+V 2>(v) - f(v)) v+1/2 - 2x[F(v^. 1/2 ) - F(v v )] 


2 \ i ' +1 ' 2 n "[ g Cv+1/2 ) - f ]. 


(6.16b) 


By convexity, the average of the two terms on the left of (6.16a) and (6.16b), 
dominates the entropy value U(»), computed at the average of the 
corresponding two Godunov solvers; invoking Theorem 6.4 we arrive at 


U[v v (t+k)] < U[v v (t)] 

- X4 [f(vv 2 ) + u '(vv 2 > ' [s v " 1/2 ) (''v-V 2 ) ■ £(v S-V 2 )] ] 

G 

+ X / V U"(g <v “ V 2)-f) + X / V /2 U M (g (v+1/ 2>- f), (6.17) 

G v 

vv 2 

and the result (6.14) follows, noting that S^ V± ^( v ^ ± l/ 2 ) = h v ± l/ 2 ’ see 
(6.7b). 

In the same way as the last lemma followed from Theorem 6.4, Theorem 6.5 


gives 
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LEMMA 6.8. 

Consider the TVD scheme (1.7a), (6.11) in the neighborhood of a non- 
critical value where Av^ y^ • Av v+ I/2 = Ttien for entr °py pairs 

(U, F) we have 


»[V t+ k>] S - kfVVz" F v- 1 / 2 ) 

V 

+ X / U"(v)[g^ V ^(v; s + = s = ~) - f(v)]dv 

G A 

W 2 

y G 

+ X / U”(v)[g^ v+ ^ 2 ^(v; s + = s = y) “ f(v)]dv. (6.18) 

v A 


Remark 6.9 . The numerical entropy flux, F^ + y^ , in (6.18) is determined by 


Godunov's intermediate value v , 


’ v+ l/^ * It coincides with the one given in 

(6.14b) and (6.14c): indeed, v G 1. = v G , 1 . (6 = V? ) is located in between 

- v+ V2 v*- 7 2 

v* + l/^ ( s+ = * *) and v v+ I /2 ( s = 1 *)• The proof is omitted. 

The only difference in Lemma 6.8 is in the slopes, s*, involved in the 

last two integrals on the right of (6.18): at a non-critical value, they are 

twice then otherwise allowed by Lemma 6.7; see (6.14a). Making use of |s^| 

to distinguish between these two cases, see (4.3a), the corresponding two 


estimates can be unified as follows. 
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THEOREM 6.10. 

Consider the difference scheme (1.7a) satisfying the TVD condition 


1 s q -v 2 S 


1 + k,l 


(6.19) 


Then, for all entropy pairs (U,F) we have 


U[v v (t+k)] < U[v v (t)] - X(F v+1/2 - F v _ 1/2 ) 


+ X f U"(v) 

vv 2 

V V+ V 2 

+ X f U"(v) 




f(v) dv 


s+ ‘ 8 ’' —ir^) ‘ £<v) 


dv, 

( 6 . 20 ) 


with a consistent numerical entropy flux 

W 2 = f (VV 2 ) " U '(V/ 2 )[VV f (vv 2 )]* (6 * 20b) 

For later purposes, we shall prefer a slightly different version of inequality 
(6.20), where the last two integrals on the right are shifted over the I y+ l/^ 


interval. We state 
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THEOREM 6.11. 

Consider the difference scheme (1.7a) satisfying the TVD condition 




‘ V v±V 2 1 


Then for all entropy pairs (U,F) we have 


"['»<»«] i "h v C‘>] - >(V>/f F v-V 2 ) 


v+1 

■ X / U"(v)| 


(v+1 h)l + 1 + K 1 - 1 + K+l , 

g 2 yi % - n — * s ,.^.- ^—r f < v > 


2X » v+1 


4 ' 


(6.22a) 


with a consistent numerical entropy flux 


F v+1/2** F <V/j) ' “'(Wz^VVz' f (vv 2 )) 

(V+ V 2 ) ( v . Hjtt l i )- f(T) 


v+1 

- X / U"(v) 

G 

VV 2 


(6.22b) 
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Remark. Theorem 6.11 finally leads us to the modified flux function ( 1) we 
were looking for; its restriction to the I y+ \/^ -interval is given, according 
to ( 6 . 22 a), by 




v+ V 2 


= g (v+ V 2 ) 


(v) 


g ( v + v 2 , ( v; .♦ . ^ . (6.23) 


Proof. We first note the identity 


U ^ 

/ /2 U"(v)g (v+1/2 ^(v; sj, s“)dv + / U"(v)g (v+1/2 ^v; s*, s 2 )dv 

Vv VV 2 

= / U"(v)g (v+1/2 ^(v; Sj, s 2 )dv. 


(6.24) 


Indeed, of the two integrals on the left, the first depends only on s^ the 
second only on s 2> Adding and subtracting 


X / V+1 U"(v)[ g <v+ 1 / 2 )(vj .♦ ■ 


1 + Ivi 1 


G 

vv 2 


f(v)|d 


to the right-hand side of (6.20a) and using (6.24), the theorem then follows. 


( 1 ) To be precise, g(v) is double—valued in the critical neighborhood, 
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Remark 6.12. The two integrals appearing on the right of (6.22) are 
responsible for the unnaturally shifted estimates w.r.t. the v^-grid 
value. If, in the previous proof, we instead add and subtract 



to the right-hand side of (6.20a) and use (6.24), we end up with a reversed 
shifted form 


U[v v (t+k)] ^ U[v v (t)] - A(F v+ i /z - F v _ 1/2 ) 

♦»/” d"c»)[s (v ' i * ) (v; - ^ 4 ^ - •; ■ ^ 41 ^) - «»>]*». 

V V-1 ' -* 


(6.25a) 


where 

W f (vv 2 ) ■ u '(V/ 2 H h v + i/ 2 - f (vv 2 )] 



(6.25b) 


Averaging of (6.22) and (6.25) results in the symmetric entropy estimate of 


this type. 
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U[v v (t+ k) ] ^ U[v v (t)] - x(F v+1/z - F v _i /2 ) 

V v V V+1 

+ X J U"(v)|g(v)j I ^ - f(v)|dv + X / U"(v)|g(v)j I ^ 


7. Necessity and Sufficiency of E-Fluxes 

Theorem 6.11 shows that a desired cell entropy inequality, (5, 
provided the last integral on the right of (6.22a) is negative 



I s J 




The requirement of negativity for all convex entropy functions 

case, is equivalent to the requirement that g(v)| T will 

v+ V2 

above f(v), depending on whether v^ v v+ ^ or ^ V v > 

account the piecewise linear form of g(v)| T in (6.23) and 

V+V2 

requirement amounts to the following inequalities: 


(1) W tVv 2 - f(v) ] 2 “• vel v+v 2 


^ ^ It is sufficient to consider 

+ + 1+ Kl 
W S --2X-> and W S 


lying in between 


f (v)jdv. 

(6.26) 


,2a), follows 

) dv < 0. 

J (7.1) 

U, in this 
lie below or 
taking into 
(6.4), this 
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(ii) (a) x Max 

|vvj < < 


f(v) - g 


1 + 1 8 J 


e v 1 + IS I C V+ V 2 ^ Av v+ V 2 ^ 


(b) 


X Max 

I'"'' u+ l I i 


| f(v) ~ Vll , 1 + l»vtll 


v+1 1 = v+1 


G v+1 1 + | s . I v+ V 2 V v+ V 2 


l AV v+l/o I* 


Numerical fluxes satisfying condition (i) are called E fluxes after Osher 
[18]. Equivalently, the corresponding schemes are characterized as exactly 
those having more numerical viscosity than that of Godunov's scheme [26] . In 
terms of the modified viscosity coefficient Q^ + y in (3.3), the first 
inequality, (i) reads, consult (2.5c), 


< + i / 2 s w <E) 

Observe that two viscosity coefficients are comparable in the sense advocated 
in [26], provided that the same modified flux is used for the corresponding 
schemes; in the above E-condition, the original flux, f^, was used in both 
cases, see Remark 3.1. Regarding the CFL condition in the second inequality 
(ii), it can be shifted to the more symmetric form in the spirit of Remark 


6 . 12 , 
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X Max 
Iv-v I < e* 
v = v 


f(v) - j 


1 + l-J 


± 2C v±vJ A VVU 

e v =■ i+iri— 


(CFL) 


With a similar shifting of the TVD condition (6.21), we arrive at 


THEOREM 7.1 (Tadmor [26, Theorem 5.1]). 

Consider an E-type TVD scheme (1.7a) such that 


A + I* I i-lvd) 

-J- 2 - •- 


and assume the following CFL condition holds 


I f (v) - g 


|v-v v | < e* I V ‘ V v 

± _ 2 C V* V? ^ Av y± V? ^ 

1 + |s | 


1 + l s J 


(7.3) 


Then, for all entropy pairs (U, F) the following cell entropy inequality is 
satisfied 

(7.4) 


U[v (t+k)] i U[v (t)] - x(F 


v+ V2 v- V2 * 


Remark. We have preferred to center the TVD condition (7.2) and CFL condition 
(7.3) around the v^-grid value. Therefore, no specific reference was given 
as to the consistent numerical entropy flux, F^ + in (7.4): it may be 
chosen, throughout the computational grid, by (6.22b), (6.25b) or any convex 
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combination of the two. 

Proof . The only part which requires clarification is the validity of the TVD 
constraint, see (6.21), 


X 


tg v+ v? 

iv v+V2 


- %+ V2 * 


(7.5) 


Indeed, (3.1) and (2.5) yield 


V V, 


vv 2 


^ l ivv 2 i + 


A v V2 


w 2 | 


Vv 2 + 


V V? 


vv 2 


(7.6) 


since g /Av , 1 , is positive by (3.7c), we find, as in (4.9b) 
v Vi V 2 


Ag 


‘v+ V? 
,v v+V 2 


Max(|g v |, |g v+1 |) 

S -——1- s 


|4v v+l/J 


‘W 2 


(7.7) 


Augmented with the E condition, (7.6) and (7.7) yield (7.5) 


x ASv+1 ^ 

Av v+ V 2 


<= Cy 2 


g V+V9 


vv 2 


< 1 


2v+ V 2 


vv 2 


^v+ V 2 * 


Remark 7.2. It is the CFL-like condition (7.3) which throws some further 
light on the first-order accuracy limitation encountered with difference 
schemes satisfying all cell-wise entropy inequalities (5.2a). Indeed, in case 
Av v + 1/ 2 ^ 0 (similarly Av v < 0), then necessarily + JL will 
lie above (similarly below) f^ in contrast to the requirement induced by 
(7.1); specifically, as v approaches v^, the CFL-like ratio on the left of 
(7.3) 
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X Max 

v-v I < 
v = 

will diverge unless the correction g^ vanishes. In this case, the modified 
flux g^ is reduced to the differential one f^, and we are back in the 
standard first-order accurate E-schemes, subject to the familiar 


f(v) - f(v ) 


(7.8) 


Max . 

-I ^ 


f(v) - f(v ) 


(7.9) 


In particular, according to (7.9), the usual half CFL number is sufficient in 
this case , in agreement with [26] 


X Max If"(v)I 

W-v v l < |Sv v±1/2 l 


(7.10) 


Of course, (7.9) offers a more delicate alternative. Note that unlike the TVD 
constraint (1.8), some further information regarding the flux behavior in 
between the grid values , is always required for a cell entropy inequality to 
hold. 

Next we show that the E condition is, in a sense described below, 
necessary for a cell entropy inequality, thus complementing the sufficiency of 
Theorem 7.1. To this end, we subtract v^(t) from both sides of (1.1) and 
(1.3), divide by k = At and take the limit At + 0, obtaining the 
semidiscrete approximation 


dv 

_v 

dt 


Ax ^v+ V 2 *V- V 2 ^ O’ 


(7.11a) 


here, the time derivative is understood in the distribution sense, and the 
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numerical flux is still allowed to depend on the mesh size, Ax, 


w 2 - h iwi''"’V £,x * °> <7 - llb> 

The corresponding cell entropy inequality will read 

§[%<*>' + s(VvVv 2 ) - 0! <7,12a> 

here, (U, F) is any consistent numerical entropy pair 

F(w,w,»«»,w; f,X + 0, Ax) = F(w) = J W IT f". (7.12b) 


In stating our next theorem, the terminology of separated data will be used: 

we shall say that the data are separated if all relevant grid values involved 

in the scheme's stencil, are separated w.r.t., say, the separating interval 

I.,i, ; that is, we either have 
J + V2 


j-k 


j+1 


j+k* 


k = l,...,p. (7.13a) 


or 


j-k 


j+1 


j+k* 


k = l,...,p. (7.13b) 


A monotone profile is, of course, the canonical example of separated data: 
any two consecutive grid values may serve as a separating interval in this 
case. 

Equipped with this terminology, we may now turn to 
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THEOREM 7.3. ( E-flux Is necessary for cell entropy Inequality) . 

Consider the semi-discrete approximation (7.11) and assume the cell 
entropy inequality (7.12) holds for all consistent numerical entropy pairs 
(U, F). 

Then the numerical flux , h, is an E-flux in any separating interval . 

Proof . Introducing the consistent numerical entropy flux 

■Ev+Va ' F<v v+i ) - "'’Vl’lVVz " (7 - 14a) 

then the following equality [18, Section 3] holds 
an i l Vy+1 

3T [v v (t)1 + aJ[W 2 "4-V 2 ] = u " (w) [VV 2 " f(w) ] dw * (7 * l4b) 

v 

Subtracting (7.12a) from (7.14b) we find, after multiplication by Ax 

Vi 

/ «w)ldw + (H^- H v _ 1/z ) < Oi (7.15) v 

V 

here. Hi, stands for the difference F . i. - F . i, , which is consistent 

v+ V 2 v+ V 2 -a)+ V 2 ’ 

with zero 

H(w,w,•••,w; f,X + 0, Ax) = F(w) - F(w) = 0. (7.16a) 

Consider a (2p+l)-grid-values stencil separated by, say, the Ij+1/ 2 interval; 
we shall extend it to all grid points, defining 

V j+k E V j+P’ V j-k 5 V j- P ’ 


k = p+1, p+2. 


(7.16b) 
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Let v be any value in between Vj and v j+i and apply (7.15)^ with 

Kruz'kov's choice for entropy function: U(w) = |w - v|: for all but the 
v = j index, the first term in the left of (7.15> v vanishes and we end up 
with 

VV 2 = H j+3/2 ^ •** ^ H » " °« v = j+1, j+2,• •• , (7.17a) 

H j-V 2 = H j-3/2 ^ ••• ^ H -„ " °» v = j-1, j-2, • • • . (7.17b) 

Here, the equalities = 0 on the right of (7.17) follow from the 

consistency (7.16a) and the definition of our extension in (7.16b). 
Complemented with (7.15)j , which in this case reads 

>1 

0 = >/ 5 < w - v >[ h j+ l/ 2 " f (w)]dw + (H j+ i/2 - 1/2 ) 

V j 

= s j+1/2 [h j+1/2 - f(v)] + (H j+1/2 - Hj . 1/2 ); (7.18) 

we find, for all v in Ij + l / 2 » 


s j + \ t V v 2 ‘ f<v) 1 * V v 2 ' V >/ 2 * °- 

That is, h, is an E-flux over the !j+1/ 2 interval as asserted. 

Next we have to confront the limitation of E schemes being at most 
first-order accurate [18]: indeed, according to the E-condition, their 
modified numerical viscosity is bounded away from that of the Lax-Wendroff 
scheme, hence (3.4) fails. Moreover, being an E-flux in separating intervals 
alone, still meets the limitation of first-order accuracy. The proof follows 
along the lines of [18], 
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8. Second-Order Resolution and a Cell Entropy Inequality 

In light of our discussion in the last section, we are seeking cell 

entropy inequalities for SOR-TVD schemes, only for special entropy pairs: 

2 

specifically, we shall consider the pair U(v) = , F(v) = f v wf'(w)dw. In 

the genuinely nonlinear case where f is, say, strictly convex, all (convex) 
entropy inequalities (and therefore convergence to the unique solution) then 
follows in the limit; see [14,5]. 

We start with 3-point E-type upwind schemes: upwinding simply means that 

we difference in the "streamwise" direction. In particular, considering the 

characteristic direction, we require either C + L i, or C - , to vanish 

v+ V 2 v+ 72 

depending on whether f'(v) is positive or negative throughout the I y+ \j^ “ 
interval; in either of these two nonsonic cases, the sum of the incremental 
coefficients i,_ equals their difference in absolute value, i.e.. 


vv 2 
W f) 


Af. 


v+ V? 

‘W 2 


f'(v)| T * o. 

V+ V2 


(8.1a) 


This is, of course, nothing but a restatement of the well-known fact that away 
from sonic values, all upwind schemes coincide with Roe-Murman's (2.5d). 

The SOR-TVD schemes to be considered are then constructed according to the 
piecewise constant viscosity modification of Theorem 4.7, where the viscosity 
coefficient used is given by 


^v+ V2 


VV 2 


VV 2 


(f) + 


A V V2 

AV v+ V 2 


(8.1b) 


Here g is the flux correction, see (4.21b), 
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Vv 2 J' a W 2 i 


with second-order resolution condition requiring 


xV) v+1/ 3 X 2 (a u+1 ) 2 + 0<|Av|) vtl < Qv+1/ (f), 


with a first-order perturbation 0(|Av|) v+ y , which is yet to be precisely 
determined. Unless otherwise stated, the following CFL-like condition will be 
assumed 

w 2 (f) ^ 1/3 * (8 - le) 

In particular, the above use of the piecewise viscosity modification from 
Theorem 4.7 is justified, since (4.15) holds. Furthermore, the essential 
estimate (4.9b) 




I 4 'V l/ 2 I ■ 2 ^ V ‘ 


yields the TVD requirement (6.21) 


48 v±l/, 




Ag v± V? 

4v v±V 2 


1 a v± v 2 1 


4 V 1/2 


. n . Ag v± V? _ n SOR , 3 , 1 

= ^v± V 2 Av u %± V 2 = 2 %± V 2 ^ = 2 1 
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so that Theorem 6.11 applies.According to this result, the desired cell 
entropy inequality holds for the above quadratic entropy provided 


! <' rf V 2 >|v; s* - * + 2 | Svi , 


is nonpositive. 

We first study the critical case, stating 


ilKsilj 


- f(v) 


dv 


( 8 . 2 ) 


Lemma 8.1. Consider the critical case |s^ | = | s v+i I = 0 ( that is, 
both v y and v v+ i are critical grid values), and assume the following CFL 
condition holds 


X Max |f'(v)| < 1/2. (8.3) 

|w v IS|4Vl/2l 


We then have 


V v+1 

/ 

V 


g ( * + 1 / 2 ) (v; 



Vi 



dv < 0. (8.4) 


Proof . We appeal to the conservative form of the modified scheme, quoted in 
(4.22): by assumption s y and therefore vanish, and hence 

Vl/2 "itVl ■ Vl/2lVlU J Since s v+l/2 agrees with the sign of 
g v+ ^, we conclude that the numerical flux correction vanisties in this 


In fact, twice the CFL limit (8.1e) can be used away from critical 
neighborhoods where [s^j • js^-jJ j £ 0. 
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case. Thus, we are left with the original E-flux we have started with. By 
Theorem 7.1 therefore, all entropy weighted integrals in (7.1) and in 
particular (8.4) are non-positive, provided the CFL condition (7.3) holds: 
the latter follows in view of (8.3) noting that, the modified flux corrections 
g^ and g v+1 vanish - consult Remark 7.2. 

Next we turn to the noncritical case |s^| + | s v+ jJ > Is calculating the 
area below the piecewise linear "double- wing" form of 


(v+i/2)^ v . s + = 1 + Kl a - = 1 + Kn[\ 


2X ’ v+1 


see figure (6.4), we obtain 


r (v+ V 2 ) / 


/ + 1+ l s vl - 1+ l s v + ll\ 

^ V ’ s v 2X * s v+l “ 2 J~ 


f(v) 


lV v+ V 2 **v+ V 2 + 2X 


-nrr[ c v + i/ 2 l 


+ 1 + is v+1 i 


v+l 


( Av v+ l/ 2 ) - / f (v) dv; 


(8.5a) 


inserting the value of h^ + from (1.7b) and (3.1), 

f + f j 1 j 

h v+ 1/2 = 2 + 2X + ®v+l ) - 2X Q v+ V 2 AV v+ V 2 ’ 


we find 
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r w 2 


g (' l+1/ 2 ) ( v . .+ 


1 + |s I 


1 + Is |\ 1 

-ir^)-H d 


‘ V v+ V 2 2 


zti - / V+1 f (v) dv 




Kf" ^ v 2 + TW (c - v 2 > 2 + (v v 2 


(8.5b) 


In the next proposition, we estimate from above the last brackets on the right 
of (8.5b). The somewhat technical proof is postponed to the end of this 
section. 


Proposition 8.2. The following estimate holds 


% * s v+1 


r v+ y 2 - Q v+ V 2 + 1 + |s„| ( C V 2 1 * 1 + I Vll (0v+ V * ] 


[Qv+i/ 2 <£) “ l2(a \tV 2 ]- 


( 8 . 6 ) 


In view of (8.5b) and Proposition 8.2, we can now complement Lemma 8.1 with 
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Lemma 8.3 . Consider the noncritical case |s | + |s + | >_ 1 ( that is. 


r (v+ V 2 


>/ ♦ 1 + i\i - 1 + J. 

\ vi S v - 2X • Vl -2X- ) ' £(v) J d 


Kti*) 2 

2X 

2\ 

K + v 2 ) 2 

» et . i/ 

v+ 1/2 

stands for 

. -interval 

h - 



l%l + I-*, I 

7 E V 1/2 *-2- IV 1/2 ' 


Vl*J 

(8.7a) 


ET V+ i /2 = Av v+ 1 /2 |- - V — y - -| " / f < V > dV * (8 * 7b > 


To ensure the integral in the left is nonpositive, hence to guarantee a 

cell entropy inequality in this noncritical case, we would therefore like to 
2 2 

choose X (a ) v+ which dominates the first two terms inside the right 
brackets of (8.7a). 


We claim that such a choice admissible by the second-order resolution 
condition (8.Id)—is indeed available in all nonsonic intervals: according to 


the upwinding property (8.1a) 
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VV 


l 2 lv‘ftl 2 


(8.8a) 


and by a standard error estimate [4, Section 2 ] 
-.3 


ET 


v+ V 2 




12 




5 e V 1 *' 


(8.8b) 


It is here that we use the first-order flexibility in defining X 2 (a 2 > v+1 ^ 
previously noticed in Theorem 4.2. We choose 


v+ V2 


%+ V2 


(f) + X 


3 


Max |f"(v)|. 
I v+ V 2 


Unifying the critical and noncritical cases, we arrive at 


xV> v+ V 2 - q; + i /2 «> I*V1 4 1 T «“* l f " (v >l - (8 ’ 9) 

v+ V2 

In order to comply with the inequality on the right of (8.Id), we set in 
(8.1c) 

g v = T- Min [ Q v± l / 2 (f 5 " x 2 (& 2 ) v± V 2 * 1 AV v± V 2 1 * (8 ' 10a) 

In terms of the numerical flux modification 


B v+ 1/2 ■ 7 [l v + Vl ' V 1/2 IS V + 1 - *,11* 


(8.10b) 


the resulting scheme then reads: 
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v v (t+k) = v v (t) - X(£ v+1 - f^,) + J [ 4 (VV2 <£,4 VV2 <t))J 

- ( K v + V 2 - K „-^) <8 - 10c> 

and we summarize what we have shown in 

THEOREM 8.4. Let Q y+ (f) be an upwind E-type viscosity coefficient, 
and assume the following CFL condition holds 


X • Max|f'(v)| < 1/3. (8.11) 

Then the difference scheme (8.9), (8.10) satisfies 

(i) Total variation diminishing. 

(ii) Second Order resolution. 

(iii) The cell entropy inequality :^ 




v+ V 2 F v+ V 2 ^ ’ 


Remark 8.5 . As argued in Theorem 4.2, the scheme has formal second order 
accuracy in all intervals such that and g^ +1 do not vanish. According 
to (8.10a) therefore, all nonsonic noncritical intervals such that 


(1) 


The numerical entropy flux is given in (6.22b) 
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Q v+1/2 «>-^ 2 ) v+1/2 - X | a v+1/2 | - X 2 | a v+1/2 | 2 

I Av , 1+ 1 /J 

- X - 5 -“- Max | f"(v) | £ 0, 

V V 2 

are included; in view of the CFL condition (8.11), the first difference on the 
right exceeds 2x|a v+ ]^ |/3, and the last inequality is valid provided that 

}l Av v+ l/ 2 l' Max ,f " (v)| = IVl/2 1, (8,l2) 

VVz 

In other words, scheme (8.9), (8.10) is second order accurate at all 
noncritical neighborhoods satisfying (8.12) . Observe that increasing 
X 2 (a 2 )^ + i^ results in a similar increase of the modified viscosity 
coefficient, 2 V+ 1/ 2 in (3.4a). 

Remark 8 . 6 . The additional viscosity added is directly related to the amount 
2 

of entropy loss across shock discontinuities which is precisely 

|ET v+ |. This should indicate the possible generalizations to other 

entropy pairs; moreover, it seems to imply that second order resolution is the 
maximum possible for cell entropy satisfying schemes, the entropy loss being 
cubic [14], 

As a special case, let us consider now the genuinely nonlinear model 
where f is, say, strictly convex. We distinguish the possibilities 

(i) The shock case. Av^ + ^ < 0. Since the trapezoidal error in 

( 8 . 8 b) as negative in this case, the admissible choice 

X 2 (a 2 )^ + = X 2 |a^ + | 2 = Q 2 + l^Cf)> will make the integral (8.7a) non¬ 

positive, as required. 
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(ii) The rarefaction case. Av . > 0. As expected, this is the more 
- v+ / 2 

delicate case because the sign of the trapezoidal error is now reversed to be 
positive, which will be compensated for by additional dissipation as before. 
To summarize, we choose 


^V) v+1/2 - Q^ 1/2 (f> (*V V J + Max [f"(v)J (8.13) 

V+ V 2 


in which case we have 


THEOREM 8.5. ( The convergence of SOR-TVD schemes) . Assume f is strictly 
convex, Q^ + (f) is an upwind E-type viscosity and the CFL condition (8.11) 

holds. Then the difference scheme (8.13) (8.10) satisfies 
(i) Total variation diminishing ; 

( ii) Second Order resolution where : 



( iii) A consistent quadratic cell entropy inequality ; 
and, as a consequence of (i) and (iii), 

(iv) Convergence to the unique physically relevant solution . 

Remark 8.6 . The above mentioned quadratic cell entropy inequality is due to 
the related inequality of Theorem 6.11. If (the nonshifted form of) Theorem 
6.10 is used instead, the same cell entropy inequality follows, this time with 
added dissipation compensating the rectangular rule error rather than the 
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trapezoidal one. Details are omitted. Recently, Osher [19] has considered 
SOR-TVD semidiscrete approximations to the genuinely nonlinear model (1.2). 
Convergence was then guaranteed with more dissipation added at shocks, rather 
than at the expected rarefactions. Indeed, the estimates used in [19] 
correspond to that quoted in our Theorem 6.10, and the dissipation unnaturally 
yet necessarily was added at shocks can be therefore attributed to the 
reversed sign of the rectangular rule error. We also note that in the cases 
studied in [19], [20] and in the last two theorems, formal second order 
accuracy had to be given up at strong jumps, which were measured w.r.t. the 
amount of convexity, Max [f"(v)]. 

I v+ V 2 

Remark 8.7 . One can obtain similar results by starting with 3-point TVD 
schemes which are not necessarily of upwind type. For example, choosing the 
modified Lax-Friedrichs schemes where Q v+ (f) = 1/2, see [26], one can 
repeat the above arguments and conclude similar convergence results. Details 


We close this section with 


Proof of proposition 8.2. The flux correction g in (8.1c) satisfies 


S tM'W” - x (>: V>/ 2 b 


- x 2 (a 2 ) * + vJ- (8 - 15> 
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We recall the definition of the viscosity coefficient in (8.16) 


v+ V2 v+ V2 


(f) 


<*y, 

4 vv 2 


Finally, we examine the incremental coefficients 
(3.1) we have 


v+ V2 


in (1.9) 


C.i /2 --y[VV 2 < f) 


4 %+ V? 

4 W 2 


+ As 1 . 

v+ *72 


A V V? 

A V V 2 


Consequently, the following equality holds 


1 ^ fc + 1 2 

1 + |s | ^v+ V 2 j 


l« I 


4^fQ v+ l /2 (f)-Xa v+1/2 ] 2 


Indeed, if Is I = 1, then both sides vanish; if otherwise s 
’ 1 v 1 V 

vanish, then 


4 ®v+ \ 

/ A ®v+ V? \ ®V+1 

*v+l 

4 V v 2 

\ 4 w 2 / 4 w 2 

4 V v 2 


vanishes as well, so that in view of (8.17a) equality holds 
Similar arguments yield 


1 ~ I 
1 + K+i 


v+lj (r - -\2 

-T ( v+1/2 3 


l S u+ll 


Q v+ V2 ^ 


' Xa v+ V 2 - 


(8.16) 

: utilizing 

. (8.17a) 

. (8.17b) 

id therefore 

in (8.17b). 

2 . (8.17c) 


In view of the TVD constraint (1.8) we also have 
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v+ V2 ^ v+ V2 2 *• \rt- V2 


2 / Ag v+ l fr \ 
\ iV v + l 


Adding (8.17b), (8.17c), and (8.l7d) we arrive at 


1 + i%i * c v+ v 2 5 + 1 + i% + ii " Cy+ v 2 ^ 


UV" + s? 


♦IfkfcJ + 


i % i - ivii)< + v 2 ‘ £ >- (8 - 18) 


The inequalities (8.15), (8.16) together with (8.18) amount to the desired 


Cv 2 (£) - l2(a V* + [- 1 + 2 Vv 2 (£> + 1st 


where in view of (8.14)^ together with the CFL condition (8.1e), the last 
brackets on the right are indeed nonpositive 


1 + 2 Vv 2 (£)+ srrrr -< 1+ ?vv 2 (£)<0; 


+ 2 Q v + l* <£) 


+ - V - V 2 -1 
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As before, the last brackets on the right are nonpositive, while the sum of 
the preceeding three terms does not exceed the asserted value 

-1V V 2 <£) + ! v 2 <£) + 1 K* v 2 (£ > - v 2 ] 

S I [Cv 2 (f) ' x 2 (a \+V 2 ] ! 

finally, although not specifically referred to, g^ = g^ + j =0 in the 
critical case, |s^| = | s v+ j| = 0, yielding 

- Vv 2 (£) + 3 Cv 2 (£) s °- 
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